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Abstract 

The atmosphere of the Sun is highly structured and dynamic in nature. From the photosphere and chromosphere 
into the transition region and the corona plasma-/? changes from above to below one, i.e. while in the lower atmosphere 
the energy density of the plasma dominates, in the upper atmosphere the magnetic field plays the governing role — one 
might speak of a "magnetic transition" . Therefore the dynamics of the overshooting convection in the photosphere, 
the granulation, is shuffling the magnetic field around in the photosphere. This leads not only to a (re-)structuring 
of the magnetic field in the upper atmosphere, but induces also the dynamic reaction of the coronal plasma e.g. due 
to reconnection events. Therefore the (complex) structure and the interaction of various magnetic patches is crucial 
to understand the structure, dynamics and heating of coronal plasma as well as its acceleration into the solar wind. 

The present article will emphasize the need for three-dimensional modeling accounting for the complexity of the 
solar atmosphere to understand these processes. Some advances on 3D modeling of the upper solar atmosphere in 
magnetically closed as well as open regions will be presented together with diagnostic tools to compare these models to 
observations. This highlights the recent success of these models which in many respects closely match the observations. 
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1. Introduction 

The solar atmosphere extends from the photo- 
sphere and chromosphere through the transition 
region into the corona. In the photosphere and lower 
chromosphere, where the Fraunhofer absorption 
lines are formed, the plasma is usually dominating 
the magnetic field, which is frozen-in. Models of the 
upper part of the convection zone and the photo- 
sphere have to account not only for the interaction 
of the plasma and the magnetic field (magneto- 
convection), but also for the r adiati ve transfer (e.g. 
iMihalas and Weibel Mihalad . Il984l) . Higher up m 
the atmosphere the plasma becomes optically thin 



* Accepted for publication in 

Advances in Space Research, 21.03.2007 

Email address: peterakis.unl-freiburg.de (H. Peter). 
Preprint submitted to Elsevier 



and the radiative transfer reduces to a radiative 
loss function, however, other processes such as heat 
conduction become of importance. 

What is common for the description of (almost) 
all interesting phenomena on the Sun is the interac- 
tion of the magnetic field and the plasma. For many 
models one assumes the magnetic field only to pro- 
vide a flow channel, e.g. in penumbral filaments or 
coronal loops. Despite the great success of many of 
these one-dimensional models one has to account 
properly for the complex magnetic structure to un- 
derstand the true nature of the solar atmosphere, 
and this implies to advance to more complex 3D 
models. 

Traditionally the 3D magnetic structure of the 
outer atmosphere has been described by extrap- 
olations of the magnetic field, from simple po- 
tential field models to very complex and compu- 
tationally intensive force-free extrapolations (e.g. 
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Schriiver et al. Such models are of vital im- 

portance to understand the large scale structure, 
but are also very helpful when, e.g., investigating 
how the magnetic field is rearranging itself. Dur- 
i ng the dynamic phases of the reconnection process 
(Buchner, 2006h . however, the field certainly be- 
comes non-force-free, and thus one also needs a 
model combining the magnetic field and the plasma, 
i.e., a MHD model. 

Due to the advancement of computer power such 
3D MHD models of structures of the solar upper 
atmosphere became available during recent years. 
Of course, these models cannot describe individual 
structures, as e.g. a coronal loop, in such great de- 
tail as ID models, but they can properly account 
for the interaction and coexistence of different struc- 
tures, which is of great importance when investigat- 
ing a highly structured object such as the solar at- 
mosphere. 

Also on global scales the interaction of mag- 
netic fields and the plasma is of gre at importance 



For ex ample, in their global models iLionello et al 



()2005l) investigate the large scale evolution of the 
(outer) corona and its connectivity to heh osphere 



(see also the review of iLinker et al.l . 120051 ) . They 



show that always different footpoints are connected 
to field lines going into the heliosphere, which is due 
to on-going reconnection. The dynamic reconfigu- 
ration processes cannot be properly handled by a 
sequence of magnetic field extrapolations. Opening 
up of loops, reconnection of open regions to form 
loops or interchange reconnection between open 
and closed regions does happen all the time. Such 
scenarios have been sketched in the past, of course, 
but only now global MHD models could show that 
they indeed operate. 

The present paper concentrates on the upper at- 
mosphere of the Sun, i.e. effects within the photo- 
sphere and chromosphere as well as on global scales 
will not be discussed. The aim of this paper is to dis- 
cuss coronal models for some structures, namely the 
closed field structures such as (moderately) active 
regions, and open structures like coronal holes. 

For illustrative purposes the structure of the low 
coronals sketched in Fig. 1. It shows the co-existence 
of closed magnetic structures of various sizes: small 
loops with lengths below 5 Mm connecting mag- 
netic patches within the chromospheric network, 
low- lying loops with length below 20 Mm span- 
ning across network cells and finally larger loops 
connecting larger magnetic patches over large dis- 
tances. The base regions of these large loops might 



be funnel-type, just like the magnetically open coro- 
nal funnels. The latter should be most prominent in 
coronal hole, but also might be present in the quiet 
predominantly magnetically closed corona. 

The paper will start with an overview of how to 
actually measure coronal magnetic fields (Sect. 2), 
even though we know only little about this currently. 
In the following Sect. 3 it will be outlined how ob- 
servations (and models) require to advance from ID 
to 3D models and why one really needs 3D mod- 
els (Sect. 4). The main part of the paper. Sect. 5, 
will discuss 3D coronal box models and the vacuum 
ultraviolet (VUV) spectra synthesized from them, 
which allows a detailed comparison to observations. 
Open structures, especially the origin and accelera- 
tion of the fast wind in coronal holes will be the sub- 
ject of Sect. 6, before Sect. 7 concludes the paper. 



2. Measurements of coronal magnetic fields 

Before discussing modeling of the solar atmo- 
sphere including the magnetic field, the attention 
should be drawn on how to actually observe the 
magnetic field in the corona. There are various 
techniques providing the coronal magnetic field. 

Direct coronagraphic measurements of the Stokes 
vector e.g. in the infrared lines of Fe XIII at 10747 A 
and 10798 A applying the Zeeman effect show mag- 
netic fiel ds in bright loo ps being of the order of 
20 Gauss ( Lin et al. . 2000[ ). Such measurements can 
be applied only for regions seen above the limb, of 
course, and are most valuable for the large-scale 
structure of the corona. Smaller areas, especially low 
in the corona, are not accessible through this tech- 
nique, because they are occulted by the design of 
the coronagraph. 

Radio observations allow to investigate the coro- 
nal magnetic fi eld also on th e disk, e.g. through the 
Zeeman effect (IWhitel . [2005l) . Radio interferometers 
such as VLA provide a remarkable spatial resolution 
of 12 arcsec, and experiments for the future FASR fa- 
cility showed that the magnetic field st r ength could 
be retrieved in a reliable way ()Whitel . [2005h . The 
problems here are the detailed interpretation that 
the emission at a given frequency does originate from 
a complex volume in the corona, and that these ob- 
servations are of value mainly for active regions. 

Another valuable tool are infrared observations in 
the He I (10830 A) triplet using Zeeman and Hanle 
effect giving magnetic fields for emerging flux re- 
gions (jLagg et al.lJ2004h . However, He I is not a coro- 
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Fig. 1. Sketch of the structure of the low corona of the Sun. This cartoon stretches some 60 Mm in the horizontal and 40 Mm 
in the vertical direction. It indicates the small loops within the network magnetic patches, low- lyin g loop s cross ing network 
cells as well as larger loops rooted in stronger magnetic patches and open coronal funnels. Following iPeted 1I2OOII) . 



nal line and can thus only indicate magnetic fields 
in just newly emerging flux regions, which still carry 
cool plasma. 

Extrapolation techniques based on photospheric 
and/or chromospheric measurements are also a very 
valuable tool to investigate coronal magnetic fields. 
As a rule of thumb one might trust these extrap- 
olation on scales typically comparable to or larger 
than a fraction of the chromospheric network, i.e., 
some 5-10 Mm. On smaller scales still the expan- 
sion of the magnetic field into the chromosphere, 
i.e. the equilibrium of magnetic and gas pressure, 
can be expected to dominate. However, whenever 
interesting events occur, such extrapolations, no 
matter how fancy, can not be trusted completely, 
as then the assumption of a force free state is most 
likely violated. While in such events the overall 
field structure on the larger scales mentioned above 
might not change dramatically, the force-free vi- 
olation might be more severe on smaller scales in 
a transient manner. This will be certainly true on 
scale s on which explosi ve events do occur (i.e. some 
Mm; llnnes et al.l . Il997l ) as well as on scales not yet 
resolvable, e.g. in the case of nanoflares. Thus such 
extrapolations are important and helpful to under- 
stand the overall structure and long-term evolution, 
but real measurements of the coronal magnetic 
field with a spatial resolution of 1 Mm and better 
are desperately needed in order to pin-point the 
relevant processes during phases when the heating 



mechanism is showing itself through a dynamic 
non-force-free event, as in flares, explosive events or 
nanoflares. While the above mainly concerns closed 
field regions, the same is also true for magnetically 
open region as found in coronal holes. F ollowing the 
idea of e.g. lAxford and McKenzid (fl997n magnetic 
reconnection of open field with bipolar network flux 
is the energy source for the heating of the open 
corona and the wind acceleration. It would be of 
high interest to detect the non-force-free magnetic 
reconfiguration predicted by this furnace model 
through the direct measurement of the magnetic 
field. 

The ideal tool to investigate the magnetic field in 
the corona is to use lines which are formed in the 
transition region and corona, i.e. to use the same 
VUV emission which is currently widely used to in- 
vestigate the corona with imagers such as TRACE 
or EIT or spectrometers such as SUMER or CDS. 
Of course, because of the complicated spatial struc- 
ture of the corona and the high variability, the in- 
terpretation of the polarized light in the VUV will 
not be possible in a simple and straightforward way, 
and is thus confronting us with the same problems as 
radio observations. The coronal models, which be- 
come more and more elaborate and realistic, will be 
of pivotal importance for a proper interpretation of 
VUV and radio spectro-polarimetric data. 

The importance of UV and VUV spectro- 
polarimetry from space has been emphasized by 
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Truiillo Bueno et al. who suggest to use the 

Lyman series of hydrogen for an investigation of the 
Hanle effect. Another approach is to measure the 
Stokes vector in strong Hnes such as the C IV dou- 
blet at 1548 A and 1550 A, which is for med in the 



transi tion region. This is suggested by iWest et al 



(|2005l ). who are currently building an instrument 
for a series of rocket flights, the first probably in 
spring 2007 for a test of the system (Solar Ultravio- 
let Magnetograph Investigation, SUMI) . Because of 
the limited count rate during the rocket flight, they 
expect to get a signal only in Stokes V, which would 
be in itself a great step forward, if they succeed. 

In the future we have to aim at a combination 
of spectro-polarimetry in the radio and the VUV. 
In the radio such techniques already exist, but they 
have to be refined, also with respect to spatial resolu- 
tion. For VUV spectro-polarimetry we will certainly 
have to wait some time until space instrumentation 
will be available, but this is the instrument to aim 
at when investigating the magnetic structure of the 
corona. 



3. Prom ID to 3D: the closed corona 

The most conspicuous basic building block of 
the corona is the coronal loop. Observations show 
that more or less semi-circular shaped loops arc the 
most prominent structure in the corona, and they 
are found in active regions, after flares, in quies- 
cent regions in the network, etc. Therefore most ID 
models for the solar corona dealt with loop models, 
one of the most wel l-known being the RTV models 



( Rosner et all . 119781 ) 



Even though it is not clear what the "micro- 
structure" of an observed loop would be, i.e. if it 
consists of many individual strands or if it i s a bun- 
dle of parallel field lines ( Klimchulj . l2006l ). many 
loop models exist solving the mass, momentum 
and energy balance along a loop-shaped ID struc- 
ture, some now with adaptive mesh refinement and 
even a self-consisten t treatment of the ion i zation 
and radiation (e. g. iBradshaw and Masonl . 12003 



Miiller et al. . 2003f ). Such loop models allow a very 
detailed description of the thermal and dynamic 
properties within the loop (or strand) modeling 
complex non-linea r processes, e.g. catastrophic 
cooling (IMiiller et a l.. 2004), which is also found in 
observations I Schriiven . i200ll ). ID models for open 
structures will be discussed in Sect. 6.1. 

The main disadvantage of these ID loop models 



is that they (usually) are not able to treat the heat- 
ing process in a physical wavPI Since the early RTV 
work almost all ID models assume a parameterized 
form of the distribution of coronal heating, e.g. ex- 
ponentially decaying with height, i.e. not really in- 
corporating a physical heating mechanism into the 
model. However some models made an attempt to 
include a distribution of heat into the ID models as 
found from other studies. 

In the mid 1970ies the Sky lab observations 
showed that the chromosphe ric network is expand- 
ing to higher temperatures ( Reeved . Il976f) . which 
lead to assume a funnel-type structure of th e tran- 
sition region and low corona (jGabriel , Il976h . How- 
ever, it became clear later that these funnels alone 
are not a good representation of the corona, even 
though they are an important ingredient (cf. Fig. 1). 
For example in the magnetically closed corona the 
emission measure EM — Jy dV, i.e. a measure 
of the emissivity of the corona at a given tempera- 
ture, could not be reproduced by the funnel models. 
Because of the ineffective heat conduction at low 
temperatures the ID models give an extremely thin 
transition region, i.e. a very steep temperature gra- 
dient, and thus in the model there is not enough 
material at low temperatures to account for the 
observed increase of the emission me a.sure towards 
lower temperatures below some 10^ K. lDowdv et al 



( 1986f ) proposed that loops at various temperatures 
could account for that increase and through this 
introduced a complicated magnetic structure of 
the low corona, sometimes also called "magnetic 
junkyard" (cf. Fig. 1). Later it was proposed that 
signatures in the spec tra of VUV emission lines 
support this scenario (|Peterj . l2000f) . Furthermore 
there is ample observational evidence from VUV 
line emission maps acquired with SUMER that low- 
lying loops with len gths below 10-20 Mm do exist 
(e.g. Feldman et al.L 2003). However, all these sce- 
narios are not really based on a physical model, but 
are merely trying to find a plausible way to draw a 
sketch of the magnetic structure of the low corona. 

In open field regions such as coronal holes it 
seems problematic that the above scenario for the 
(predominantly) closed field regions holds. There 
other proce s ses m ight be more important. Recently 
lEsser et al.l ( 20051 ) have shown that the (fast) solar 
wind outflow of the plasma along funnel type struc- 



^ This is not true for the I D corona and wind mo dels of 
e.e. lMarsch and Tul lll997bl ') or iTu and MarschI lll997l) as dis- 
cussed in Sect. 6.1. 
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tures might well result in significant violation of the 
ionization equilibrium by lifting a sizable fraction of 
neutral material higher up into hotter regions of the 
funnel. This can then produce an increased emis- 
sion in Ly-Q; which agrees well with observation. 
However, it still has to be explored to which extent 
this process will predict the observed emission from 
the transition region, i.e. of lines formed at a couple 
of 100 000 K. 

While large-scale structures such as a (polar) 
coronal hole or an active region complex seem to 
be rather stable for days and weeks, on smaller 
scales of below 10 Mm the corona is highly dynamic. 
When analyzing the spatial magnetic structure on 
these smaller scales this variability is adding to the 
complexity of the problem. For example, when mod- 
eling an explosive event, it is usually assumed that 
the "background" magnetic structure remains the 
same, while only a small perturbation is shuffling 
around the magnetic field le a ding t o a reconnection 
event (e.g. Innes and Tothl . 1999^ ). When investi- 
gating observations with imagers or spectrographs, 
one sees a very strong variability in time indicating 
that the "natural" state of the the transition region 
and corona o n scales below 10 Mm is a dynamic one 
(|lnnesl . l2004l ). 



4. Why 3D models? 

One could argue that a composite model consist- 
ing of a large number of loops and open structures 
would be sufhcient to describe the solar corona. Fur- 
thermore one could argue to neglect the magneti- 
cally open structures in these models as they provide 
only a minor contribution to the emission in X-rays 
or the far VUV, i.e. in the wavelength bands accessi- 
ble to Yohkoh, or the EIT or TRACE coronal chan- 
nels. As the corona and transition region are built up 
by (loop- like) magnetic structures which are rooted 
in the solar photosphere, one could use an extrapo- 
lation of the magnetic field to identify the loops and 
then investigate the resulting structure. 

For certain proble ms this is certain l y a go od ap- 
proach. For example ISchriiver et al.l (|2004[ ) did a 
global extrapolation of the observed photospheric 
magnetic field to define the loops and used a sim- 
ple scaling law for the flux of mechanical energy Fh 
into the loop, Fh oc f. Here B is the mag- 

netic field at the base of the structure, L its half 
length, / ^ exp[— (i?/500 G)^] a factor accounting 
for reduced heating in sunspots, and /? and A are 



free parameters. Each loop was then described in a 
static ID model. The authors then calculated how 
this constructed multi-loop corona would appear in 
an X-ray observation (Yohkoh SXT) and compared 
it to the real observation. By varying /? and A they 
found that a heating scaling with (3 ~ 1.0±0.5 and 
A w — 0.7±0.3 would give the best fit. Even the best 
fit is not really representing the real corona, which is 
not surprising as it is based on a very simple model. 
Nevertheless it can be used to estimate how the en- 
ergy input into the corona is roughly scaling with the 
magnetic structures. Most important, this kind of 
model opens the possibility for more realistic global 
models of stellar corona, which are not accessible to 
direct observations. 

For other purposes one-dimensional models, or 
such multi-loo p models would n o t be sufficient. 
For example iSchriiver and Title! (|2003l ) showed 
that depending on the distribution of internetwork 
magnetic flux only some 50% of the coronal mag- 
netic field above the quiet Sun does come from 
network magnetic patches. This implies that the 
simple idea of a coronal loop being rooted in a 
single (or even a few) patches of strong magnetic 
flux i n the photosphere is not ap plicable. Further- 



more iJendersie and Peter! ([20061) could show that 
based on current instrumentation we are not even 
able to determine where the coronal field is really 
rooted, on a scale of a fraction of a super-granule 
(say 5 Mm) . For only small differences in the dis- 
tribution of weak internetwork flux elements one 
gets radically different connections from the photo- 
sphere into the corona. Therefore when describing 
the corona in more detail on a scale corresponding 
to a fraction of a single active region or the chromo- 
spheric network, one has to account for the complex 
magnetic connectivity, because a simple picture 
of a coronal structure rooted in a single (or few) 
strong magnetic patch(es) does not really apply. On 
top of that the problem is even more complicated, 
because the distribution of magnetic flux in the 
photosphere changes constantly on time scales from 
minutes (granules, ~l Mm) to days (super-granules, 
w20 Mm), which is shorter than the lifetime of many 
loops. Loop arcades, e.g. as seen by TRACE, overall 
last longer (many days) than the nano-flaring pro- 
cesses associated with them, so they appear to be 
rigid, or solid. However, the magnetic connectivity 
within the loop structures might constantly change 
while the loop is sitting seemingly still, i.e. the 
many strands of the loop could constantly change 
their identity. Currently it is not clear if, and if yes 
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to what extent, this affects the appearance of the 
loop as a whole. Therefore we might ask ourselves 
the question what sense the concept of a coronal 
loop as a rigid magnetic loop really makes. 

In consequence, to describe the structure of the 
transition region and corona, one has to account for 
the changing "boundary conditions" in the photo- 
sphere and include a proper treatment of the 3D 
structure of the magnetic field. 



5. 3D MHD coronal box models 

One of the most important ingredients for a coro- 
nal model is a proper energy equation. Because of 
the strong flows observed in the active regions as 
well as quiet corona, the enthalpy flux can play a 
major role. And only if the energy balance, including 
furthermore the energy input, heat conduction and 
radiation are treated well enough, the density of the 
corona will be correct (as good as possible...). This 
is because basically the energy input sets the pres- 
sure of the coronal structure: for larger energy input 
more energy will have to be radiated in the transition 
region. When the heating rate increases this conse- 
quently leads to (chromospheric) evaporation, while 
the coronal temperature changes only little as the 
heat conduction acts as a thermostat through the 
strong dependency on temperature (ocT^/^). Thus 
only a proper energy equation will assure proper 
densities and consequently proper velocities within 
the computational domain. 

A severe problem is that the energy dissipation 
processes are occuring on scales much smaller than 
currently observable or resolvable with the MHD- 
type models to describe e.g. an active region complex 
or a coronal loop. Thus in terms of kinetic physics 
it is questionable is a "proper energy equation" in 
fluid terms can even be defined (e.g. Marscb, 2006). 
However, one might (or have to) hope that when 
describing the corona with a resolution of the cur- 
rent MHD models (some km in ID, some 100 km in 
3D), the energy deposition as a function of space and 
time is comparable to the result a micro-physical 
model smoothed to the MHD scales would give. Also 
the resistivity in the MHD models is basically set 
by the (spatial) resolution they can achieve. Due to 
the limitations for the magnetic Reynolds number 
this is a common property for all MHD-type mod- 
els. However MHD studies with different spatial res- 
olutions and resistivities showed that the energy de- 
position does not change significantly over a range 



of scales (for resolut ion and resistivity) access i ble to 
MHD experim ents (Gals gaard and Nordlund . 1996t 
Hendrix et al. |,[l996). Future work will have to show 
to what extent the MHD approach is a good rep- 
resentation for the real Sun, but currently it is the 
most appropriate (and only) tool to describe struc- 
tures such as active region systems or loop com- 
plexes. 

The basic test for every coronal model has to be 
whether it is matching with observations. The most 
straight-forward way for such a test is to compare 
observable quantities such a line intensities or shifts 
synthesized from the model to real observations, 
as this eliminates the many problems of inversions 
of y UV spectral data (e.g. Judge and McIntoshL 
In the following a 3D MHD coronal box for- 
ward model will be discussed which compares very 
well with observations concerning e.g. the emission 
measure, line shifts or temporal variability. 



5.1. 3D forward box model 



A 3D MHD model presented bv lGudiksen and Nordlund 
(j2005af dh includes the atmosphere from the photo- 
sphere to the lower corona in a 60x60 Mm horizon- 
tal times 37 Mm vertical box. It solves the mass, 
momentum and energy bal ance, th e latte r including 
classical heat conduction ( Spitzer . 19561 ) and opti- 
cally thin radiative losses as piecewise power laws. 
The temperature of the chromosphere is kept near 
a prescribed profile by Newtonian cooling. Initially 
the magnetic field is given by a potential field ex- 
trapolation with the lower boundary as observed by 
MDI/SOHO for an active region. Further on in the 
simulation the field evolves self-consistently and be- 
comes non-potential, of course. The system is driven 
by horizontal motions at the lower boundary. The 
flow field is constructed using a Voronoi-tessellation 
technique ( Okabe et al. . 1992[ ) and reproduces the 
typic al pattern of the granular motions of the Sun 
(Schrijver et al.l . ri997h . By this also the power spec- 
tra of the velocity and the vorticity are reproduced. 
The spacing of the (non-uniform) grid goes down to 
150 km and is as good it can be given the size of the 
structure modeled here (60 Mm). Of course all the 
shortcomings of MHD models as discussed above 
apply here. For the upper boundary condition it 
is assumed that the magnetic field above the com- 
putational domain is potential, which, of course, 
prevents the development of larger open magnetic 
structures (and was not the aim of that study) . This 
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has to be kept in mind when discussing the Doppler 
shifts of coronal Unes (Sect. 5.2) and the apparent 
lack of hot blueshifted material, i.e. a wind. 

This procedu re leads to a heating of the corona 
just in the way IParked (|l972 t) suggested. The field 
line braiding gives high coronal temperatures of a 
million K, and a system of hot loops is forming, con- 
necting the magnetic concentrations of the active 
region. Furthermore the system reaches some sort 
of quasi-stationary state, with large fluctuations in 
time and space, or in other words, the system is 
quite dynamic. The detaile d discussion of the MHD 
results can be found in Gudiksen and Nordlund 



can 

(|2005al lbh. 

Based on the above M HD model r e sults f or the 
density and temperature Peter et al. ( 2004 . 2006f l 



calculated the emissivity at each grid point under 
the assumption of ionization equilibrium using the 
atomic data base CHIANTI. Assigning a spectral 
line profile at each grid point with a line width cor- 
responding to the thermal width, one can integrate 
the spectra from the 3D box along a line of sight, 
e.g. the vertical. This results in 2D maps of spec- 
tral profiles, which can be analyzed in the same way 
as observations, e.g. through maps of line intensi- 
ties or shifts. For a detailed discussion of the synthe- 
sis of the emission lines from the transition region 
and corona and the assu mptions and limitations see 
Peter et all (|2004l2006l ). 

Figure 2 shows in the middle panel the vertical 
magnetic field at the bottom of the computational 
box, i.e. in the photosphere. To the left and right are 
maps in line intensity and line shift as derived from 
the spectra computed from the MHD model. They 
represent a view from the top onto the box, corre- 
sponding to an observation at disk center. As found 
with observations the spatial structures in the tran- 
sition region line CiV (1548 A) at 10^ K are much 
smaller and finer than in the corona seen in MgX 
(625 A) formed at lO^K. 



5.2. Diagnostics of the 3D model: VUV spectra 

In the following a short overview will be given 
how VUV spectra synthesized from a 3D coronal 
box model can be used to test the model. It should 
be stressed that the model results presented in the 
following show the best overall match to the obser- 
vations published so far for the differential emission 
measure, Doppler shifts or rms fluctuations. 
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Fig. 3. Differential emission measure (DEM) from a 3D MHD 
box model compared to observations. The solid line shows 
the fit from the DEM inversion based on the lines displayed as 
bars, which have been synthesized from the MHD model (in- 
tegrated over the whole box). The thick dashed line is based 
on a DEM inversion using observed quiet Sun disk center 
line radiances observed with SUME R llWilhelm et al.L|f998ah 
scaled by a factor of two. Following iPeter et al. I iI20o4k 



Diffential emission measure 



For each time step iPeter et al.l (|2004l ) used the 



radiances of a number of VUV emission lines syn- 
thesized from a MHD model to perform a differ- 
ential emission measure (DEM) analysis, i.e. they 
used the maps of the spectral profiles ("synthetic 
observations" ) as described in the previous section 
and integrated the line profiles over the whole map. 
Thus the line radiances represent the total emission 
of the given VUV lines integrated over the whole 
computational domain. As with real observations 
one can then use an inversion procedure to de- 
rive the differential emission measure (DEM), e.g. 
with the help of the procedures provided with the 
CHIANTI atomic data package (jPere et all . Il997 : 
Young et al.l . l2003l) F1 Fig. 3 shows the resulting 
DEM curve as a solid line for a single time step 
derived from the intensities of a given set of lines, 
which characterizes the (average) small active re- 
gion simulated in the MHD model. For comparison 
the thick dashed line shows the inversion of actually 
observed intensities of the same lines. The match of 
the model with the observations is rem arkable, es- 



pecial ly, the 3D MHD model used in the [Peter et al 



(|2004h study l|Gudiksen and Nordlundl . l2002h 



is 



^ As with any DEM inversion, many implicit assumptions 
apply, e.g. ionization equilibrium, constant abundances, con- 
stant pressure atmosphere. 
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Fig. 2. Maps in line intensity and Doppler shift in the transition region (left two panels) and the corona (right two panels) as 
derived from a snapshot of a 3D MHD coronal model. In these maps the computational box is viewed from above, i.e. this 
represents the situation when observing near disk center. The middle panel show s the vertical magnetic field at the bottom of 
the computational domain, i.e. in the photosphere. Following IPeter et al.l l|2006h . 
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Fig. 4. Average Doppler shifts of the synthesized spectra 
when viewing the computational box from straight above for 
two different time steps (bars and rectangles) . The height of 
the bars indicate the standard deviation of the Doppler shifts 
for the respective time step. The thick dashed lin e show s 
the trend found in observations jPeter and Judgj. Il999l) . 
Following IPeter et al.l H2004V 



reproducing the observed trend of the emission 
measure below ^10^ K, while previous (ID and 
2D) models failed to predict this increase of DEM 
towards lower temperatures (cf. Sect. 3). 

In the above model the increase towards the low 
transition region is caused by numerous low-lying 
(intermittent) co ol structures. This co nfirms the sce- 
nario outhned bv lDowdv et al. ( 19861 ). in which a hi- 
erarchy of large and small, hot and cool loops exist in 
the transition region and corona as already discussed 
in Sect. 3. Despite the variability of the low corona, 
the emission measure (averaged over the computa- 
tional domain, viz. the active region) changes only 
little with time. 



Average Doppler shifts 

Since the discovery of the systematic tr ansition 
region line rcdshifts bv lDoschek et al.l (|l976i ) , the ex- 
planation of these persistent Doppler shifts has been 
a challenge for modelers. The discovery that the 
lin es form e d in th e (lo w) corona show a n e t blue shift 
by IPeteil (|l999l ) and IPeter and Judgel (|l999f) us- 
ing SUMER/SOHO data added a new quahty to 
this challenge. The trend of the net Doppler shift 
at quiet Sun disk center is shown i n Fig. 4 as a 



thick dashed line ( d ata co r npiled from .Brckke et al 



19971: IChae et al.l , 
iTeriacaet al.Lll999l ). 



19981 IPeter and Jud ge. ,199S : 



In order to derive the average Doppler shift cor- 
responding to disk center observations one can in- 
vestigate the Doppler maps of the synthetic spectra 
with a vertical line-of-sight, i.e. when looking at the 
computational box from straight above (cf. Fig. 2 



The a verage line shifts as computed by IPeter et al 



(|2004l ) for two different time steps 7min apart are 
plotted in Fig. 4 as bars and rectangles (the heights 
of the bars and rectangles indicating the scatter of 
the Doppler shifts)!^ 

The line shifts are caused by the dynamic response 
of the atmosphere to the energy input due to the 
braiding of magnetic field lines. The flows are partly 
induced through energy deposition in and subse- 
quent expansion of the corona, partly from evapo- 
ration of chromospheric material due to heating at 
low heights and the heat flux from the corona into 
the chromosphere. 



Actually, the average values for the shifts of the Doppler 
maps are quite similar to the Doppler shifts of the average 
spectra. 
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Fig. 5. RMS fluctuations in line intensity and shift from the 
spectra derived from the MHD coronal model for a number 
of VUV lines as a function of line formation temperature. 
T he thick dashed line sh ows the tre nd found in observ ations 
bv lBrkovic et al ] l|2003h . Following IPeter et al 

Throughout the transition region (<5-10^ K) the 
match of the observed Doppler shifts with those 
of the spectra synthesized from the model is quite 
good. It should be stressed here that no fine-tuning 
was applied, but these Doppler shifts follow nat- 
urally from the driving of the corona through the 
footpoint motions of the magnetic field lines. In the 
low corona (>5-10^K) the synthesized spec tra do 
not s how the blueshifts as they are observed (j Peter , 
I999I ) , but this problem might be overcome (at least 
partly) by a more appropriate handling of the upper 
boundary condition. 



Temporal variability 

Based on the synthesized maps of line intensities 
and shifts generated from a line of sight integra- 
tion along the vertical direction, the rms fluctua- 
tions at each spatial pixe l within the r espec tive map 
have been evaluated by IPeter et all (|2005l) . Figure 
5 shows their results for a number of VUV lines as 
a function of line formation temperature. The dia- 




Mg X (625 A) 



relative intensity log(//</>) 



relative intensity log(//</>) 



Fig. 6. 2D histograms of plasma-/3 as a function of line 
intensity (a; CIV, b: MgX). The dashed lines show the 
median value of plasma— /3 in a given interval of intensity. 
Please no te that the 2D histo grams are plotted on a log-scale. 

Following iPeter et al. I Hooi). 

monds show the (spatially) averaged values for the 
rms fluctuation for each line and the bars represent 
the scatter of the rms values (standard deviation) . 

This is exa c tly th e same procedure as used by 
Brkovic et all (|2003[ ) ;o reduce their observational 
data from CDS and SUMER on SOHO. The trends 
they found with their observations is over-plotted as 
thick dashed lines in Fig. 5. It has to be emphasized 
that the spatial resolution of the maps in line in- 
tensity and shift derived from the MHD model has 
been reduced to match the spatial resolution of the 
SOHO instruments, as this has an effect on the ab- 
solute values of the rms fluctuations (the trend with 
temperature, however, does not depend on the spa- 
tial resolution). 

Together with the good match between the aver- 
age Doppler shifts and emission measure from ob- 
servations and the spectra computed from the 3D 
MHD coronal box model presented above, this is yet 
another strong indication that the underlying 3D 
coronal model is a good description for the coronal 
heating mechanism. 

5.3. The magnetic structure 

To investigate the interaction of the magnetic field 
and the plasma one can discuss plasma-/3, i.e. the 
ratio of the gas pressure to the magnetic pressure, 
(3 — 2iip/B^ . Fig. 6 shows the relation of (3 to the 
synthesized transition region (a; C IV) and coronal 
(b; Mg x) emission within t he computational box as 
done bv lPeter et all (|2006f) . The plots show the 2D 
histograms of the plasma beta as a function of (nor- 
malized) intensity of the respective line, i.e. the color 
coding represents the fraction of the volume where 
one finds the respective combination of plasma-/? 
and normalized intensity. The median relations be- 
tween plasma-/? and the normalized intensities are 



9 



over-plotted as dashed lines. The transition region 
and corona are mostly in a low-/3 state, however with 
noticeable exceptions in the coronal part! 

It is clear that in the whole transition region the 
assumption of a low-/3 plasma is very good (Fig. 6a). 
This is different, however, in the corona when relat- 
ing (5 to the Mg X emission (Fig. 6b) . A large fraction 
of the not-low-beta plasma in the corona is at low 
densities, i.e. at log //(/) < 0, where the magnetic 
field is very weak, too. As these are regions of low 
emission, one might be inclined to disregard this. 
However, even in the volume with a relative intensity 
larger than the median intensity (i.e. log //(/} > 
in Fig. 6b) contributing ^90% to the total coronal 
Mg X emission, so me ~5% of this bri ght material has 



values of /3 > 1 ( Peter et al. . 2006f ). These are re- 



gions of rather high densities but low magnetic field, 
and they can be found e.g. in regions along magnetic 
neutral lines across which the direction of the mag- 
netic field changes and thus strong currents heat the 
plasma and finally cause the density to increase. 

In conclusion, in significant patches of the corona, 
even above active regions, (3 is in fact not much less 
than unity. 

This again emphasizes that for "good" coronal 
model it is not enough to only extrapolate the mag- 
netic field and then solve a lot of ID loop-like prob- 
lems along each field line. Instead, one has to account 
for the interaction of the plasma and the magnetic 
field in a more complex 3D model. 



6. 3D MHD models for open structures 

So far the discussion concentrated on magneti- 
cally closed structures, such as active regions or the 
chromospheric network. In magnetically open re- 
gions the plasma is not trapped by the magnetic 
field, but can escape along magnetic field lines. This 
is the basic reason why coronal holes, i.e. the source 
region of the fast solar wind, appear dark. Because 
a large part of the energy goes into acceleration of 
the wind, there is less energy to heat the plasma. 
Consequently the temperature is a bit lo wer than in 



the qu iet Sun, typically below 10 K (Wil helm et al 



1998bl ). The effect of the reduced heating rate is 
much stronger on the pressure, which results in a 
lower density in the corona above the holes — we 
see less emission from the coronal holes. 



6.1. ID scenarios and models for the fast wind 

As for the quiet and active Sun structures, which 
are dominated by loops, the first models for the open 
corona, which is the same as for the fast wind, were 
ID models stretching from the base of the corona 
into interplanetary space, with an area expansion 
factor accounting for the (super) radial expansion of 
the magnetic field. 

As with the loops these models used a parame- 
terized form of the energy input, which was typi- 
cally concentrated close to the Sun. As mentioned 
earlier, in order to properly understand the mass 
loss, especially to describe the mass loss rate self- 
consistently set by the energy input, one has to in- 
clude a transition re gion to account fo r the heat 
flux b ack to the Sun (|Hammeij . Il982al lbl: IWithbroeL 
1988h . Only then one can describe properly the in- 
terplay between coronal heating, densities, temper- 
atures and solar wind acc eleration, as done first by 
Hansteen and Leer ( 19951 ). However, as for the ID 
loop models (cf. Sect. 3), these models did not in- 
clude a physi cal heating mechanism. 

In contrast. lAxford and McKenzie ( 1997 ) put for- 
ward an idea that small-scale reconnection in the 
magnetic concentrations of the chromospheric net- 
work (i.e. between open and closed field lines) leads 
to high-frequency Alfven waves which propagate up- 
wards into the corona, a concept also called mag- 
netic furnace. These waves can then resonantly in- 
teract with the protons and heavy ions in the open 
corona and drive the w ind a s has been descr i bed b y 
Marsch and T\3 (|l997bh and iTu and MarschI (|l997l ). 
Furthermore one can also account for the expansion 
of the magnetic field directly above the magnetic 
concentr ations in the chromosphere forming large 



funne ls (jMarsch and Tul . Il997al : iHackenberg et al 
20001 ). 



These ID models including a physical mechanism 
for the heating and acceleration of the plasma de- 
scribe a situation where the material is continuously 
and constantly accelerated f rom the chro mosphere 
out into the solar wind (e.g. Holzei . lioOSl ). By this 
they implicitly assume that the magnetic configu- 
ration from the upper chromosphere into the solar 
wind is rather stable! This would strongly contra- 
dict the scenario outlined for the closed field regions 
in the previous section, where the magnetic field in 
the photosphere is constantly shuffled, resulting in 
a continuous change also higher up. 
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6.2. Observing the source of the fast wind 

It has been known for a while that the fast wind 
is originating from coronal holes, and observations 
of blueshifts in coronal lines within coronal holes 



confi rmed this (jRottman et al.L 119821 : lOrrall et al 
19831 ). Even the transition region lines of CiV 
showed a higher fraction of blueshifte d regions in 



coronal holes tha n in the quiet Sun (jPere et al 
19891 : IPeteij . 119991) . 



The strongest bl ueshifts are found in the dark- 
est coronal regions ( Wilhelm et al. . 200Ci l and they 
to be concentrated a t inte rsections of net- 
19991) . 



seem 
work 



lanes ( Hassler et al 



comparison 



of the Doppler shift patterns to the photospheric 
magnetic field showed that these strong blueshifts, 
interpreted as the base of the outflow, coincide 
with strong mag netic concentrations in the network 
(|Xia et al.l . 120041 ). These data seemed to support the 
wind scenario described above, namely that funnel- 
type structures emerging from strong network mag- 
netic patches define the flow channels of the wind 
reaching out into the interplanetary medium. 

One can investigate the magnetic connectivity 
from the photosphere to the corona in coronal holes 
and compare this to the Doppler shift signals of 
transition region and (low) coronal lines. This shows 
that the blueshifts in coronal lines coincide with 
the magnetic funnels, but that at those locations 
no systematic blueshifts can be found in transition 
regio n lines within the instrumental uncertain- 
ties! TuetaL , . [iooi). Based on this observation 



Tu et al. (|2005^ argued that the plasma of the solar 



wind outflow has to be injected somewhere between 
the formation heights of the transition region and 
coronal lines they inspected (CiV and NeVIIl), 
i.e. somewhere around 10 Mm. If this is correct, it 
changes dramatically our view of the formation of 
the wind at its very base, because it would challenge 
the old paradigm that there is a continuous outflow 
of the wind from the chromosphere into the helio- 
sphere. However, because of the higher densities in 
the transition region the expected outflow speed 
of the wind should be rather low (wlkm/s) at the 
level where one can expect C IV to form. This is just 
at the resolution limit of current instrumentation, 
e.g. SUMER/SOHO can determine Doppler shift s 
down to some l-2km/s (jPeter and Judge! . Il999l ). 
Therefore it remains unclear, whether there is a 
continuous outflow from the chromosphere, with a 
small not yet detectable velocity at CiV heights, or 




Fig. 7. Initial setup of the magn etic field in the 3D box 
model for an open field region bv lBiichner and Nikutowskil 
1 I2OO5I) . The bottom boundary of the magnetic field is taken 
from observations and is driven by whirl flows, which brings 
toget her the open and closed regions le ading to reconnection, 
From lBiichner and Nikutowsk 



jg ions le a 
J 1I2OO5I) . 



if the velocity in C IV is really zero and the wind is 
injected bet ween the l evels of CiV and NeVIII as 
suggested bv iTu et al.l (|2005l ). 



6.3. Accounting for the magnetic structure: 3D 
models for the onset of the fast wind 



To investigate the scenario proposed by iTu et al 



(l2005h a numerical expe r iment was conducted by 



Biichner and Nikutowski ( 2005f) . They start with a 



magnetic configuration given at the lower bound- 
ary from observations and then drive the system by 
whirl flows imposed on the lower boundary. Through 
this they bring together open and closed magnetic 
field structures which results in reconnection. Their 
model contains a plasma-neutral gas coupling, i.e. 
it includes effects beyond MHD. For the reconnec- 
tion they used a switch-on-resistivity, i.e. the re- 
sistivity is non-zero only where the current-carrier- 
velocity (ocj/n) is above a certain threshold. By this 
their treatment of the reconnection is more detailed 
than the 3D box models for the closed corona dis- 
cussed in Sect. 5.1, however, they used a simplified 
energy equation, not accounting for heat conduction 
or energy losse s throu gh optically thin radiation (cf. 
iBiichner etall . I2OO5I ). The magnetic configuration 
(at the initial state) is shown in Fig. 7. 

As a consequence of driving the magnetic field 
through the whirls at the bottom, the reconnection 
between the open and closed regions injects energy 
and plasma into the open funnels and through this 
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drives a wind out of the funnels. Roughly, the plasma 
below 10^ K falls down, plasma above 10^ K leaves 
the computational box through the top. The down- 
drafts are mostly concentrated above the magnetic 
concentrations at the bottom boundary. Thus this 
sir aulation suppor ts the general picture suggested 
bv iTu et aL ( 2005 ). However, further modehng with 
a proper inclusion of an energy equation and also an 
investigation of the VUV emission to be expected 
from this model is needed before a final conclusion 
can be drawn. 

At the moment we have to conclude that even in 
the very quiet coronal holes, the assumption of a sta- 
tionary magnetic field configuration simply channel- 
ing the solar wind outflow is not justified. To really 
understand the origin of the solar wind we have to 
account for the complex interaction between differ- 
ent magnetic field structures and also between the 
magnetic field and the plasma. 



7. Conclusions 

A realistic model of the (upper) solar atmosphere 
has to account for the complexity of its magnetic 
structure. Together with the constant rearrange- 
ment of the magnetic field at photospheric levels this 
leads to an ongoing energy release in the low parts 
of the transition region and corona and through 
this drives the highly dynamic upper atmosphere. 

The current box models for magnetically closed 
regions as discussed in Sect. 5 show a qualitatively 
good match to observed quantities, such as emis- 
sion measures, Doppler shifts of rms- variability. This 
presents ample evidence that the heating through 
field line braiding is a prime candidate to heat the 
corona. Of course, other means of energy input into 
the corona might also reproduce the average quan- 
tities of e.g. emission measure, Doppler shift or line 
widths as presented in this paper. It remains to be 
seen if other models based on different heating pro- 
cesses will give the same or different results, once 
they are pushed to produce observables as now done 
for the field line braiding model. In the future more 
efforts have to be undertaken not only to increase 
the spatial resolution of the simulations, but also to 
actually include a physical mechanism for the dissi- 
pation of magnetic energy during the reconnection 
events. 

The box models for magnetically open region, i.e. 
for the onset of the fast wind (Sect. 6), have not yet 
been analyzed in such detail concerning observable 



properties, but they support some recent interpre- 
tation of observational data suggesting that plasma 
is injected in funnel-type structures to form the fast 
solar wind. 

Together with the very successful global coro- 
nal models, which have only briefly been touched 
upon in this paper, now complex three-dimensional 
models become available accounting properly for 
the plasma and the magnetic field as well as their 
interaction. We might hope that the coming years 
will bring further advances of these models and 
their diagnostic capabilities. Thus we will improve 
our understanding of the structure, dynamics and 
heating of the upper atmospheres of the Sun and 
solar-like stars. 
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